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INTRODUCTION 


Our  goal  is  to  develop  expert  systems  software  for  solving  time-dependent 
partial  differential  equations.  The  software  should  allow  users  to  describe 
problems  in  a  natural  language,  have  a  convenient  geometric  description  inter¬ 
face,  and  not  require  knowledge  of  sophisticated  numerical  analysis.  The 
systems  should  be  intelligent,  efficient,  reliable,  robust,  and  able  to  solve  a 
large  class  of  problems  to  prescribed  error  tolerances. 

The  power  of  adaptive  techniques  is  that  they  are  capable  of  making  deci¬ 
sions  that  change  the  computational  environment.  This  significantly  minimizes 
the  number  of  a  priori  decisions  demanded  of  the  user  and  provides  dramatic 
savings  in  the  cost  of  the  computation.  This  capability  is  performed  by  proce¬ 
dures  that  monitor  intermediate  results  and  feed  this  data  back  to  a  control 
mechanism  that  modifies  the  solution  strategy.  Three  popular  adaptive  tech¬ 
niques  for  solving  partial  differential  equations  are  mesh  moving  or  rezoning 
(r-ref inement) ,  mesh  refinement  (h-ref inement) ,  and  order  enrichment 
(p-ref inement) .  In  r-ref inement,  the  mesh  is  moved  either  continuously  or  stat¬ 
ically  at  discrete  times  in  order  to  resolve  nonuniformities  and  reduce  errors. 
H-ref inement  involves  the  addition  or  deletion  of  computational  cells  to  the 
mesh  and  p-refinement  involves  increasing  or  decreasing  the  order  of  a  method  in 
different  portions  of  the  domain.  All  strategies  attempt  to  organize  the  com¬ 
putation  so  that  little  effort  is  expended  in  regions  where  the  solution  is 
smooth  and  a  much  greater  effort  is  devoted  to  regions  where  the  solution  is 
more  difficult  to  compute. 

The  different  refinement  strategies  are  being  combined  to  yield  remarkable 
results.  Babuska  and  Szabo  (ref  1)  showed  that  an  hp-ref inement  scheme  produced 
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an  exponential  rate  of  convergence  on  a  singular  elasticity  problem.  Arney  and 
Flaherty  (ref  2)  developed  an  hr-ref i nement  scheme  that  moved  a  'base'  coarse 
mesh  to  follow  important  dynamic  structures  of  the  solution  and  recursively 
refined  the  base  mesh  to  improve  resolution.  They  found  that  mesh  motion  was 
inexpensive  relative  to  mesh  refinement  and  reduced  dispersive  errors  associated 
with  wave  motion  but  did  not  always  accurately  follow  structures,  especially 
when  interactions  occurred,  and  could  not  dependably  satisfy  prescribed  toler¬ 
ances.  Recursive  mesh  refinement  can  satisfy  prescribed  tolerances  but  involves 
more  complicated  data  structures  and  greater  care  at  coarse-fine  mesh  interfaces 
than  r-ref i nement . 

There  are  numerous  other  variations  of  the  three  adaptive  strategies  for 
time-dependent  problems.  For  example,  temporal  refinement  can  be  done  globally 
to  produce  an  adaptive  method  of  lines  strategy  (refs  3,4)  or  locally  in  com¬ 
bination  with  the  spatial  refinement  strategy  (refs  5,6). 

Accurate  a  posteriori  error  estimation  is  essential  for  codes  that  strive 
to  satisfy  user-prescribed  error  tolerances.  Error  estimation  is  often  the  most 
expensive  part  of  an  adaptive  algorithm.  Arney  and  Flaherty  (ref  2)  calculated 
the  local  discretization  error  at  nodes  of  the  mesh  using  an  algorithm  based  on 
Richardson's  (ref  7)  extrapolation.  This  pointwise  estimate  can  then  be  used  to 
construct  several  global  measures  of  the  discretization  error.  The  advantage  of 
this  method  is  that  it  can  be  used  to  find  error  estimates  for  any  numerical 
scheme  without  explicitly  knowing  the  exact  form  of  the  error.  Details  of  this 
error  estimate  and  its  implementation  on  a  moving  mesh  are  discussed  in  Arney 
(ref  8)  and  Arney  et  al.  (ref  9). 

In  this  report,  we  apply  Arney  and  Flaherty's  (ref  2)  adaptive' mesh  moving 
and  refinement  technique  to  one-dimensional  hyperbolic  systems.  As  described 
herein,  their  approach  consists  of  moving  a  base  mesh  of  quadrilateral  cells  to 


isolate  important  spatial  structures  of  the  solution.  Refinement,  when  needed, 
is  performed  within  cells  of  coarser  meshes.  Solutions  are  generated  by  a 
MacCormack  (ref  10)  finite  difference  scheme,  and  local  error  estimates,  which 
are  used  to  control  mesh  motion  and  refinement,  are  computed  by  Richardson's 
(ref  7)  extrapolation.  Our  goal  is  to  quantify  the  relative  costs  and  benefits 
of  mesh  motion  and  local  mesh  refinement.  In  this  report,  we  present  the 
results  of  computational  experiments  performed  on  three  one-dimensional  problems 
using  several  conventional  and  adaptive  numerical  procedures.  The  results 
obtained  demonstrate  both  the  potential  and  limitations  of  the  adaptive  algo¬ 
rithm.  We  have  mixed  results  showing  that  the  effects  of  mesh  moving  can  be 
problem-dependent.  Generally,  mesh  motion  is  effective  for  following  an  iso¬ 
lated  structure,  but  much  less  so  when  structures  interact.  In  this  report,  we 
also  discuss  the  utility  of  our  methods,  the  computational  results,  and  future 
work . 

ALGORITHM 

We  consider  an  application  of  Arney  and  Flaherty's  (ref  2)  adaptive  proce¬ 
dure  to  one-dimensional  vector  systems  of  hyperbolic  conservation  laws  having 
the  form 

ut  +  fx(x,u,t)  =  0,  x  e  D,  t  >  0  (1) 

u (x , 0 )  =  Uq ( x ) ,  x  e  D  U  3D  (2) 

with  appropriate  well-posed  conditions  on  the  boundary  3D  of  a  domain  D.  Like 
them,  we  discretize  Eqs.  (1)  and  (2)  using  a  MacCormack  (ref  10)  finite  dif¬ 
ference  scheme  because  of  its  general  applicability  (ref  11).  Although  this 
scheme  suffers  a  reduction  in  order  on  a  moving  nonuniform  grid,  Our  com¬ 
putations  show  that  proper  mesh  moving  can  provide  enough  efficiency  and 
accuracy  to  compensate  for  this  order  reduction. 
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The  MacCormack  scheme  produces  spurious  oscillations  near  discontinuities 
because  it  is  a  centered  scheme  with  second  order  accuracy  on  a  uniform  mesh. 

The  use  of  artificial  viscosity  to  make  this  scheme  total  variation  diminishing 
( TVO )  makes  it  attractive  as  a  general  solver  for  problems  with  discontinuities 
and  we  use  a  model  by  Davis  (ref  12).  The  artificial  viscosity  terms  are  calcu¬ 
lated  from  the  solution  data  at  the  beginning  of  each  time  step  and  are  added  to 
the  solution  after  the  MacCormack  solution  has  been  calculated. 

Arney  and  Flaherty's  (ref  13)  mesh  moving  procedure  is  based  on  an 
intuitive  approach  that  allows  nodes  to  follow  local  nonuniformities  rather  than 
the  more  analytical  approaches  of  equidistribution  of  error  (ref  4)  or  the 
solving  of  variational  problems  to  minimize  some  given  functional  (ref  14)  which 
can  be  expensive  and  problem-dependent.  They  derive  equations  for  the  nodal 
velocities  so  t-hat  the  mesh  moves  to  follow  the  geometric  propagation  of  some 
local  nonuniformity.  This  generally  reduces  dispersive  errors  and  allows  the 
use  of  larger  time  steps  while  maintaining  accuracy  and  stability.  Important 
factors  for  m<»sh  moving  are  to  maintain  mesh  smoothness  by  controlling  adjacent 
cell  ratios,  to  keep  nodes  within  the  domain  boundaries,  and  to  move  nodes  with 
a  velocity  that  reduces  discretiza4--;on  error.  In  order  to  prevent  mesh  distor¬ 
tion  that  can  lead  to  increased  discretization  error  of  the  solver,  mesh  points 
cannot  move  independently  but  must  be  coupled  to  at  least  some  of  their  neigh¬ 
bors  . 

Some  schemes  do  this  coupling  by  attraction  and  repulsion  of  nodes  (cf. 

Rai  and  Anderson  (ref  15)).  In  these  algorithms,  the  coupling  is  done  globally, 
where  each  node  influences  the  velocity  of  all  other  nodes  in  the  mesh. 
Attempting  to  equidistribute  errors  can  lead  to  problems  where  nodes  move 
incorrectly  in  some  regions.  This  occurs,  for  example,  when  a  mesh  that  is 
following  one  structure  must  react  to  another  nonuniformity  that  arises  in 
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another  part  of  the  domain.  An  abrupt  grid  adjustment  can  be  eliminated  if  the 
influence  is  more  local  and  the  movement  algorithm  is  combined  with  a  mesh 
refinement  scheme  to  add  tne  necessary  nodes  in  the  region  of  the  new  structure. 

At  each  time  step,  the  selection  mechanism  of  Arney  and  Flaherty's  (ref  2) 
mesh  moving  algorithm  uses  the  current  node  locations  and  the  nodal  values  of  a 
mesh  movement  indicator  at  the  independent  moving  nodes  of  a  coarse  mesh  as 
feedback.  The  local  error  estimates  are  used  as  the  mesh  movement  indicators. 
Nodes  with  'significant  error'  are  grouped  into  error  clusters.  This  clustering 
separates  the  important  spatially  distinct  phenomena  of  the  solution.  As  time 
evolves,  the  clusters  can  move,  change  size,  collide,  or  separate.  At  each  time 
steo,  new  clusters  can  be  created  and  old  ones  can  vanish. 

Mesh  movement  is  then  determined  by  each  node's  relationship  to  its  nearest 
error  cluster  and  the  propagation  velocity  of  the  center  of  error  mass  of  the 
cluster.  Therefore,  the  nodal  influence  is  regional.  The  amount  of  movement  is 
determined  by  a  movement  function  which  insures  that  the  center  of  error  of  the 
cluster  moves  according  to  a  differential  equation  suggested  by  Coyle  et  al. 

(ref  16) 

r  +  Xr  =  0  ( 3 ) 

where  r(t)  is  the  position  of  the  center  of  error  mass  of  a  cluster  and  [’)■■- 
d()/dt.  Additionally,  this  movement  function  smoothes  the  mesh  motion  and  pre¬ 
vents  nodes  from  moving  outside  the  domain  boundary.  The  distance  a  node  moves 
is  reduced  near  boundaries  to  prevent  it  from  leaving  the  domain.  Nodes  on  the 
domain  boundaries  are  not  allowed  to  move. 

Arney  and  Flaherty  (ref  2)  perform  static  rezoning  whenever  computation 
with  the  current  mesh  would  be  counterproductive  or  when  the  current  mesh  suf¬ 
fers  from  poor  mesh  ratios.  There  are  sophisticated  algorithms  to  check  the 
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mesh  condition  and  to  verify  the  validity  of  the  mesh  (cf.  Babuska  and  Szabo 
(ref  1)  and  Simpson  (ref  17)).  These  algorithms  can  check  for  gaps  between 
cells  and  overlapping  cells.  Since  moving  meshes  can  only  develop  such  severe 
problems  over  time,  mesh  degradation  can  be  discovered  before  it  develops 
complete  invalidity.  This  mesh  degradation  or  i 1 1 -condi t ioning  occurs  when  the 
mesh  angles  are  severe,  the  mesh  contains  poor  mesh  ratios  or  poor  aspect 
ratios,  or  the  time  step  is  too  restrictive  because  of  the  crowding  of  nodes. 
Nodes  of  the  coarse  mesh  can  become  too  crowded  when  error  clusters  pass  through 
boundaries  or  when  two  or  more  error  clusters  converge  and  trap  nodes  between 
them.  Static  rezoning  is  performed  only  when  absolutely  necessary  due  to  the 
high  cost  in  accurately  interpolating  the  solution  from  the  existing  mesh  onto  a 
new  one.  It  was  not  done  in  any  of  the  examples  presented  in  this  report. 

Arney  and  Flaherty  (ref  2)  used  local  mesh  refinement  to  insure  that  the 
user-prescribed  error  tolerance  was  satisfied.  This  was  done  by  recursively 
introducing  finer  meshes  by  binary  refinement  of  space-time  cells  in  regions 
where  nodes  with  unacceptable  error  have  formed  clusters  (cf.  Berger  (ref  18), 
Flaherty  and  Moore  (ref  6),  and  Gropp  (ref  19)).  The  clustering  algorithm  used 
for  refinement  is  the  same  as  the  one  used  for  mesh  movement.  The  clusters  are 
buffered  so  that  high  error  nodes  are  in  the  interior  of  the  refined  region. 

The  problem  is  recursively  solved  on  these  fine  meshes  until  the  error  is  within 
the  specified  tolerance.  The  refined  subgrids  that  are  adaptively  created  by 
the  local  refinement  algorithm  overlay  the  coarser  grids.  Each  of  these 
subgrids  is  independently  defined.  Figure  1  shows  a  coarse  grid  with  portions 
overlayed  by  two  fine  grids  and  three  finer  grids. 

Arney  and  Flaherty's  (ref  5)  mesh  refinement  strategy  suggests  the  use  of  a 
tree  data  structure  for  its  description  and  implementation.  In  this  tree  struc¬ 
ture,  the  coarsest  grid  is  the  root  node  and  is  defined  as  level  0  in  the  tree. 
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The  subgrids  of  the  coarse  grid  are  its  offspring  in  the  tree  and  are  defined  as 
level  1.  A  grid  at  level  i  is  properly  nested  in  the  tree  between  its  parent  at 
level  i  -  1  and  its  offspring  (if  any)  at  level  i  +  1.  Grids  at  the  same  level 
are  given  an  arbitrary  ordering.  Oue  to  the  clustering  and  buffering  of  error 
regions,  grids  at  the  same  level  of  a  two-dimensional  problem  can  intersect  and 
overlap.  Figure  1  depicts  an  example  of  a  sequence  of  meshes  that  might  be  pro¬ 
duced  by  our  refinement  procedure  for  a  coarse  grid  refinement  step.  The  num¬ 
bers  next  to  the  grids  indicate  the  order  in  which  the  solution  is  computed  on 
each  grid. 
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Figure  1.  Typical  set  of  local  refinement  grids  for  one  coarse  time  step. 

The  numbers  indicate  the  order  in  which  the  solution  is  computed 
on  each  grid. 


7 


Such  tree  data  structures  are  commonly  used  in  adaptive  mesh  refinement 
orocedures  (cf.  Berger  and  Oliger  (ref  20)  and  Flaherty  and  Moore  (ref  6)). 
Additionally,  we  use  a  stack  to  implement  the  recursive  algorithm  (cf.  Aho  et 
al.  (ref  21)  and  Horowitz  and  Sahni  (ref  22)). 

The  solution  vectors,  error  estimates,  and  nodal  information  are  all  stored 
in  a  dynamic  storage  area  with  pointers  from  the  tree  to  this  storage  area  for 
each  mesh  in  the  tree.  For  each  grid,  we  store  its  level  in  the  tree  and  the 
number  of  nodes  it  contains.  The  dynamic  storage  area  contains  the  solution 
vector  and  the  error  estimate  at  each  node,  and  nodal  information  for  use  in  the 
solver  and  grid  interface  procedures.  Since  the  old  mesh  data  is  saved  to 
obtain  initial  data  for  newly  refined  grids  and  nodes  of  the  parent  grids  are 
updated  from  the  fine  grid  solutions,  nodal  relationships  between  meshes  are 
stored  directly  in  the  nodal  vector. 

COMPUTATIONAL  RESULTS 

We  conducted  experiments  of  Arney  and  Flaherty's  (ref  2)  adaptive  mesh 
strategy  using  three  problems  and  our  results  follow.  In  each  case,  errors  are 
measured  in  the  norm  and  the  CPU  times  are  normalized  to  unity.  All  calcula¬ 
tions  were  performed  on  an  IBM  30810  computer. 

Example  1 

Consider  the  scalar  hyperbolic  differental  equation 

ut  +  (cos  )Tt)  ux  =  0,  t  >  0,  -0.4  £  x  $  1.4  (4) 

with  initial  conditions 

I  1,  if  0.4  <  x  ^  0.6 

u(x,0)  =  j  (5) 

0,  otherwise 

and  boundary  conditions 

u ( -0 . 4 , t )  =  u( 1 . 4 , t )  =  0.  (6) 
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The  exact  solution  to  this  problem  is 


u(x,t)  =  < 


(7: 


1,  if  0.4  $  x  -  (sin  7tt)/rr  $  0.6 
0,  otherwise 

which  is  a  square  pulse  of  unit  amplitude  that  oscillates  sinusoidally  about  the 
center  of  the  domain.  Artificial  viscosity  was  added  to  eliminate  oscillations 
in  the  solution;  however,  this  resulted  in  an  attentuation  and  spreading  of  the 
square  pulse. 

Four  different  adaptive  strategies  were  used  to  solve  this  problem  for 
0  $  t  $  2.5.  The  solutions  at  several  times,  the  mesh  trajectories,  and  the 
time  step  profile  for  the  various  strategies  are  shown  in  Figures  2,  3,  4,  and 
5.  Table  I  summarizes  the  comDutational  cost  and  accuracy  of  the  four  strate¬ 
gies. 


TABLE  I.  COMPARISON  OF  THE  DIFFERENT  ADAPTIVE  STRATEGIES  FOR  EXAMPLE  1 


Adaptive  Strategy 

n 

Number  of 
Space-Time  Cells 

Normalized 
CPU  Time 

Attentuation 

1 .  Stationary 
uniform  mesh 

0.1090 

774 

1.000 

0.545 

2.  Moving  mesh 

0.0903 

1134 

1.452 

0.730 

3.  Stationary 
uniform  mesh 
with  refinement 

0.0614 

15718 

8.761 

0.969 

4.  Moving  mesh 
with  refinement 

0.0395 

16554 

10.069 

0.994 

With  a  stationary  uniform  mesh,  we  find  that  the  square  pulse  is  rapidly 
attenuated  and  diffused.  The  time  step  profile  shows  how  the  Courant  number  is 
utilized  to  maintain  maximum  step  sizes  without  loss  of  stability.  From  Figure 
3,  it  is  apparent  that  the  results  improve  when  the  mesh  is  allowed  to  move. 
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The  pulse  -is  attenuated  less  and  the  error  is  reduced,  but  more  time  steps  are 
needed  to  complete  the  computation.  The  mesh  trajectories  in  Figure  3  demon¬ 
strate  how  well  the  nodes  track  the  sauare  oulse  as  it  oscillates.  Figures  4 
and  5,  depicting  the  results  of  strategies  3  and  4,  respectively,  show 
remarkable  improvement  when  adaptive  mesh  refinement  is  used.  In  both  cases, 
the  local  error  tolerance  was  specified  as  0.001.  Errors  are  reduced  and  the 
attenuation  of  the  oulse  is  almost  negligible,  but  shape  distortion  is  still 
significant.  Notice  how  well  the  refinement  procedure  tracks  the  oulse; 
however,  the  cost  of  computation  increases  by  almost  an  order  of  magnitude. 

When  moving  and  refinement  are  combined,  the  results  are  even  more  remarkable. 
The  oulse  is  attenuated  by  a  factor  of  only  0.6  percent. 

Example  2 

Consider  the  linear  uncouoled  system 
ut  +  (cos  rrt)ux  =  0, 

t  >  0,  -0.4  S  x  $  1.4  (8) 

Vt  -  (cos  7Tt)vx  =  0, 


with  initial  conditions 

u(x,0)  =  v(x,0) 

and  boundary  conditions 


1,  if  0.4  $  x  $  0.6 

0,  otherwise 


(9) 


u ( -0 . 4 , t )  =  u ( 1 . 4 , t )  =  v(-0.4,t)  =  v ( 1 . 4  ,  t )  =  0.  (10) 

The  exact  solution  to  this  problem  is 


1, 

— J- 

— h 

O 

/A 

X 

1 

(sin 

7T  t  )  /  7T  $  0.6 

0, 

otherwi se 

1, 

if  0.4  $  x  + 

(sin 

7T  t )  /  7T  $  0.6 

0, 

otherwise . 

The  first  comoonent  u  is  the  same  as  in  Example  1,  and  the  second  component 
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v  moves  symmetrically  with  u.  Four  different  adaptive  strategies  were  used  to 
solve  this  problem  for  0  «  t  $  1.5.  Table  II  summarizes  the  computational  cost 


and  accuracy  of  the  four  strategies.  The  solutions  at  several  times,  the  mesh 
trajectories,  and  the  time  step  profile  for  mesh  strategies  3  and  4  are  shown  in 
Figures  6  and  7,  respectively. 

It  is  clear  that  mesh  moving  does  not  Drovide  the  expected  imorovement  in 
the  results  for  this  Droblem.  In  fact,  we  can  see  from  Table  II  that  each  time 
the  mesh  is  moved,  the  error  in  the  comouted  solution  increases.  This  is 
oecause  two  identical  error  regions  moving  symmetrically  about  the  center  of  the 
domain  do  not  contribute  equally  to  the  mesh  motion  due  to  asymmetries  in  their 
error  estimates.  As  a  result,  the  mesh  moves  incorrectly  and  the  solution 
deteriorates.  This,  in  turn,  leads  to  further  imbalance  of  the  error  clusters 
and  subsequently  causes  catastrophic  effects.  Comparing  Figures  6  and  7,  we  see 
how  bad  the  solution  is  attenuated  due  to  incorrect  mesh  motion.  Improper  mesh 
motion  has  also  led  to  refinement  in  some  regions  of  the  mesh  where  it  should 
not  have  been  necessary.  In  both  cases,  the  local  error  tolerance  was  SDecified 
to  be  0.005. 


TABLE  II.  COMPARISON  OF  THE  DIFFERENT  ADAPTIVE  STRATEGIES  FOR  EXAMPLE  2 


Adaptive  Strategy 

n 

Number 
Space-T ime 

of 

Cells 

Normal i zed 

CPU  Time 

i . 

Stationary  uniform 
mesh 

0.1145 

1650 

1.000 

2. 

Moving  mesh 

0.1221 

5640 

3.386 

3. 

Stationary  uniform 
mesh  with  refinement 

0.0541 

20828 

6.926 

4. 

Moving  mesh  with 
refinement 

0.0583 

48954 

18.667 
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Example  3 

Consider  the  coupled  hyperbolic  system  from  the  wave  equation 
ut  -  vx  =  0, 

t  >  0,  -0.3  ^  x  $  1.4 

vt  “  ux  = 

with  initial  conditions 


u(x,0)  = 


1,  if  0.4  $  x  $  0.6 

0,  otherwise 


1,  if  0.5  C  x  $  0.7 

v(x,0)  = 

0,  otherwise 

and  boundary  conditions  satisfying  the  exact  solution 

u(x,t)  =  (p(x+t)  +  q(x-t))/2.0 
v ( x , t )  =  ( p ( x+t )  -  q(x-t))/2.0 

where  , 


2, 

if  0.5  $  $ 

$  0.6 

pm  *  i. 

if  0.4  «  £ 

<0.5  or  0.6  <  £  $  0.7 

(18 

0, 

otherwise 

1, 

if 

0.4  $  $ 

$  0.5 

-1, 

if 

0.6  «  ? 

$  0.7 

(19 

0, 

otherwise. 

aU)  = 


Four  different  adaptive  strategies  were  used  to  solve  this  problem  for 
0  $  t  ^  0.6.  Table  III  summarizes  the  computational  cost  and  accuracy  of  the 
four  strategies.  The  solutions  at  several  times,  the  mesh  trajectories,  and  the 
time  step  profile  for  mesh  strategies  3  and  4  are  shown  in  Figures  8  and  9, 
respectively. 

Once  again,  mesh  motion  does  not  appear  to  result  in  the  desired  improve¬ 
ment  in  the  solution.  In  this  case,  there  are  two  error  regions  moving  away 
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with  unit  speed  in  opposite  directions  from  the  center  of  the  domain.  However, 
the  error  regions  are  not  identical  as  was  the  case  in  Example  2.  With  a  moving 
mesh,  the  solution  is  attenuated  and  consequently,  the  error  measure  in  the 
norm  increases. 


TABLE  III.  COMPARISON  OF  THE  DIFFERENT  ADAPTIVE  STRATEGIES  FOR  EXAMPLE  3 


Adaptive  Strategy 

m 

Number  of 
Space-Time  Cells 

Normal i zed 
CPU  Time 

1.  Stationary  uniform 

0.1141 

930 

1.000 

mesh 

2.  Moving  mesh 

0.1057 

3180 

3.454 

3.  Stationary  uniform 

mesh  with  refinement 

0.0527 

23236 

11.493 

4.  Moving  mesh  with 

refinement 

0.0552 

39234 

20.232 

CONCLUSION 

We  have  experimented  with  the  adaptive  mesh  method  of  Arney  and  Flaherty 
(ref  2).  Our  results  indicate  that  proper  mesh  moving  can  efficiently  reduce 
errors.  However,  their  mesh  moving  is  not  effective  for  problems  that  have  more 
than  one  moving  structure.  We  find  that  whenever  there  are  two  error  regions, 
the  mesh  moving  strategy  is  unable  to  make  an  accurate  decision.  This  occurs 
particularly  during  the  time  when  the  two  structures  have  not  completely 
separated  but  still  form  one  large  error  cluster.  Results  of  local  refinement 
tests  show  that  it  can  efficiently  reduce  errors.  The  most  powerful  method  was 
the  combination  of  both  mesh  moving  and  mesh  refinement.  Results  obtained  for 
Example  1  show  that  a  totally  adaptive  mesh  strategy  can  be  extremely  effective. 
The  overhead  associated  with  the  clustering  and  dynamic  data  structures  is  only 
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about  5  percent  of  the  time  needed  to  calculate  a  comparable  solution  on  a  uni¬ 
form  mesh. 

Additional  computation  is  needed  to  verify  the  generality  of  these  conclu¬ 
sions.  It  is  also  not  clear  how  many  of  the  difficulties  were  due  to 
MacCormack's  (ref  10)  finite  difference  scheme  or  Richardson's  (ref  7) 
extrapolation-based  error  estimate.  A  TVD  scheme  would  greatly  improve  perform¬ 
ance  near  discontinuities. 

We  are  re-examining  the  entire  process  in  order  to  determine  an  effective 
mesh  procedure.  Future  computations  will  be  performed  using  more  advanced  shock 
capturing  difference  schemes  (e.g.,  Engquist  and  Osher  (ref  23)'. 
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Figure  2.  Solutions  at  t  =  0,  0.57,  1.15,  1.80,  and  2.33, 
mesh  trajectories,  and  time  step  profile  for 
strategy  1  of  Example  1. 
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3.  Solutions  at  t  =  0,  0.87,  1.42,  1.88,  and  2.37, 
mesh  trajectories,  and  time  step  profile  for 
strategy  2  of  Example  1. 
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Figure  4.  Solutions  at  t  =  0,  0.57,  1.15,  1.80,  and  2.33, 
mesh  trajectories,  and  time  steD  profile  for 
strategy  3  of  Example  1. 
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gure  5.  Solutions  at  t  =  0,  0.77,  1.37,  1.91,  and  2.50, 
mesh  trajectories,  and  time  step  profile  for 
strategy  4  of  Example  1. 
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Solutions  for  u  at  t  =  0,  0.17,  0.48,  0.96,  and  1.38, 
mesh  trajectories ,  and  time  step  profile  for  strategy 
3  of  Example  2. 
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Figure  7.  Solutions  for  u  at  t  =  0,  0.76,  0.95,  1.17,  and  1.49, 
mesh  trajectories,  and  time  step  profile  for  strategy 
4  of  Example  2. 
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Figure  8.  Solutions  for  u  at  t  =  0,  0.20,  and  0.58,  mesh  trajectories, 
and  time  step  profile  for  strategy  3  of  Example  3  with  error 
tolerance  =  0.005. 
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Figure  9.  Solutions  for  u  at  t  =  0,  0.41,  and  0.60,  mesh  trajectories, 
and  time  step  profile  for  strategy  4  of  Example  3  with  error 
tolerance  =  0.005. 
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DIRECTOR 

US  ARMY  INDUSTRIAL  BASE  ENGR  ACTV 
ATTN:  AMXIB-P  1 

ROCK  ISLAND,  IL  61299-7260 

COMMANDER 

US  ARMY  TANK-AUTMV  R&D  COMMAND 
ATTN:  AMSTA-DDL  (TECH  LIB)  1 

WARREN.  MI  48397-5000 

COMMANDER 

US  MILITARY  ACADEMY  1 

ATTN:  DEPARTMENT  OF  MECHANICS 
WEST  POINT,  NY  10996-1792 

US  ARMY  MISSILE  COMMAND 
REDSTONE  SCIENTIFIC  INFO  CTR  2 

ATTN:  DOCUMENTS  SECT,  BLDG.  4484 
REDSTONE  ARSENAL,  AL  35898-5241 

COMMANDER 

US  ARMY  FGN  SCIENCE  AND  TECH  CTR 
ATTN:  DRXST-SD  1 

220  7TH  STREET,  N.E. 

CHARLOTTESVILLE,  VA  22901 

COMMANDER 

US  ARMY  LABCOM 

MATERIALS  TECHNOLOGY  LAB 

ATTN:  SLCMT-IML  (TECH  LIB)  2 

WATERTOWN,  MA  02172-0001 


NOTE:  PLEASE  NOTIFY  COMMANDER,  ARMAMENT  RESEARCH,  DEVELOPMENT,  AND  ENGINEERING 
CENTER,  US  ARMY  AMCCOM,  ATTN:  BENET  LABORATORIES,  SMCAR-CCB-TL , 
WATERVLIET,  NY  12189-4050,  OF  ANY  ADDRESS  CHANGES. 
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NO.  OF 
COPIES 

COMMANDER 

US  ARMY  LABCOM,  ISA 

ATTN:  SLCIS- IM-TL  1 

2800  POWDER  MILL  ROAD 
ADELPHI,  MD  20783-1145 

COMMANDER 

US  ARMY  RESEARCH  OFFICE 

ATTN:  CHIEF,  IPO  1 

P.O.  BOX  12211 

RESEARCH  TRIANGLE  PARK,  NC  27709-2211 
DIRECTOR 

US  NAVAL  RESEARCH  LAB 

ATTN:  MATERIALS  SCI  &  TECH  DIVISION  1 

CODE  26-27  (DOC  LIB)  1 

WASHINGTON,  D.C.  20375 


NO.  OF 
COPIES 

COMMANDER 

AIR  FORCE  ARMAMENT  LABORATORY 
ATTN:  AFATL/MN  1 

EGLIN  AFB ,  FL  32542-5434 

COMMANDER 

AIR  FORCE  ARMAMENT  LABORATORY 
ATTN:  AFATL/MNF 

EGLIN  AFB,  FL  32542-5434  1 

METALS  AND  CERAMICS  INFO  CTR 
BATTELLE  COLUMBUS  DIVISION 
505  KING  AVENUE 

COLUMBUS,  OH  43201-2693  1 


NOTE:  PLEASE  NOTIFY  COMMANDER,  ARMAMENT  RESEARCH,  DEVELOPMENT,  AND  ENGINEERING 
CENTER,  US  ARMY  AMCCOM,  ATTN:  BENET  LABORATORIES,  SMCAR-CCB-TL , 
WATERVLIET,  NY  12189-4050,  OF  ANY  ADDRESS  CHANGES. 


